Fluctuation-Dominated Phase Ordering Driven by Stochastically Evolving Surfaces 
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We study a new kind of phase ordering phenomenon in coarse-grained depth models of the hill- 
vaUey profile of fluctuating surfaces with zero overall tilt, and for hard-core particles sliding on such 
surfaces under gravity. We find that several such systems approach an ordered state with large 
scale fluctuations which make them qualitatively different from conventional phase ordered states. 
We consider surfaces in the Edwards- Wilkinson (EW), Kardar-Parisi-Zhang (KPZ) and Golubovic- 
Bruinsma-Das Sarma-Tamborenea (GBDT) universality classes. For EW and KPZ surfaces, coarse- 
grained depth models of the surface proflle exhibit coarsening to an ordered steady state in which 
the order parameter has a broad distribution even in the thermodynamic limit, the distribution of 
particle cluster sizes decays as a power-law (with an exponent 6), and the 2-point spatial correlation 
function has a cusp (with an exponent a = 1/2) at small values of the argument. The latter 
feature indicates a deviation from the Porod law which holds customarily, in coarsening with scalar 
order parameters. We present several numerical and exact analytical results for the coarsening 
process and the steady state. For linear surface models with dynamical exponent z, we show that 
a = {z — l)/2 for z < 3, a = 1 for z > 3, and there are logarithmic corrections for z — 3, 
implying a — 1/2 for the EW surface and 1 for the GBDT surface. Within the independent 
interval approximation we show that a + 9 — 2. We also study the dynamics of hard-core particles 
sliding locally downwards on these fluctuating one-dimensional surfaces and flnd that the surface 
fluctuations lead to large-scale clustering of the particles. We find a surface-fluctuation driven 
coarsening of initially randomly arranged particles; the coarsening length scale grows as ~ t^^^ . The 
scaled density-density correlation function of the sliding particles shows a cusp with exponent a ~ 
0.5, and 0.25 for the EW and KPZ surfaces. The particles on the GBDT surface show conventional 
coarsening (Porod) behavior with a ~ 1. 
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I. INTRODUCTION 



Phase ordering dynamics describes the way in which 
domains of an ordered state develop when an initially 
disordered system is placed in an environment which pro- 
motes ordering. For instance, when a simple ferromagnet 
or alloy is quenched rapidly from very high to very low 
temperatures T, domains of equilibrium low-T ordered 
phases form and grow to macroscopic sizes. A quantita- 
tive description of the ordering process is provided by the 
time development of the two-point correlation function; 
asymptotically, it is a function only of the separation 
scaled by a length which increases with time, typically as 
a power law ||l|. 

New phenomena and effects can arise when we deal 
with phase ordering in systems which are approaching 
nonequilibrium steady states. In this paper, we study a 
coupled-field nonequilibrium system in which one field 
evolves autonomously and influences the dynamics of 
the other. The system shows phase ordering of a new 
sort, whose principal characteristic is that fluctuations 
are very strong and do not damp down in the thermo- 
dynamic limit — hence the term fluctuation-dominated 
phase ordering (FDPO). 

In usual phase ordering systems such as ferromagnetic 



Ising model, if one considers a finite system and waits 
for infinite time, then the system reaches a state with 
magnetization per site very close to the two possible val- 
ues of the spontaneous magnetization, nis or — m^., with 
very infrequent transitions between the two. This is re- 
flected in a probability distribution for the order param- 
eter which is sharply peaked at these two values, with 
the width of the peaks approaching zero in the thermo- 
dynamic limit (Fig. |l](a)). By contrast, in the FDPO 
steady state, the system continually shows strong fluc- 
tuations in time without, however, losing macroscopic 
order. Accordingly, the order parameter shows strong 
variations in time, reflected eventually in a probability 
distribution which remains broad even in the thermody- 
namic limit (Fig. 0(b)). 

The physical system we study consists of an indepen- 
dently stochastically fluctuating surface of zero average 
slope, on which reside particles which tend to slide down- 
wards guided by the local slopes of the surface. Some- 
what surprisingly, a state with uniform particle density 
is unstable towards large scale clustering under the ac- 
tion of surface fluctuations. Eventually it is driven to 
a phase-ordered state with macroscopic inhomogeneities 
of the density, of the FDPO sort. Besides exhibiting a 
broad order parameter distribution, this state shows un- 
usual scaling of two-point correlation functions and clus- 
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ter distributions. It turns out that much of the physics of 
this type of ordering is also captured by a simpler model 
involving a coarse-grained characterization of the surface 
alone, and we study this as well. A brief account of some 
of our results has appeared in M . 
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FIG. 1. Schematic depiction of Prob(m) against m in 
steady state for (a) a normal phase ordering system such as a 
ferromagnet at low temperature (6) a system showing FDPO. 



In the remainder of the introduction, we first dis- 
cuss the characteristics of FDPO vis a vis normal phase- 
ordered states. We then discuss, in a qualitative way, the 
occurrence of FDPO in the surface-driven models under 
study. The layout of the rest of the paper is as follows. In 
Section |l| we define and study the coarsening and steady 
states of three different coarse-grained depth models of 
the fluctuating surfaces. In Section [II, we demonstrate 



the existence of a power-law in the cluster size distribu- 
tion, and show how it can give rise to FDPO. In Section 
IV, we discuss ordering of sliding particles on fluctuat- 
ing surfaces. In Section 0, we explore the robustness of 
FDPO with respect to changes in various rates defining 
the nonequilibrium process. Finally, in Section VI 



we 



summarize our principal results, and discuss the possible 
occurrence of FDPO in models of other physical systems. 



A. Ordered States in Equilibrium Systems 

With the aim of bringing out the features of fluctuation 
dominated phase ordering (FDPO) in nonequilibrium 
systems, let us recall some familiar facts about phase or- 
dered states in equilibrium statistical systems. We first 
discuss different characterizations of spontaneous order- 
ing, following the paper of Griffiths on the magneti- 
zation of idealized ferromagnets. We follow this with a 
discussion of fluctuations about the ordered state. 



1. Definitions of Spontaneous Order 

(a) In the absence of a conservation law, the magnetiza- 
tion m is an indicator of the ordering: 



(1) 



where L is the linear size, d is the dimension and s„ is 
spin at site n. In the thermodynamic limit, the thermal 
average of the absolute value 



mi = LinriL- 



(nonconserved) 



(2) 



with Boltzmann-Gibbs weights for conflgurations pro- 
vides an unequivocal measure of the order. This is be- 
cause in the low-temperature ordered phase, the prob- 
ability Prob^m) of occurrence of magnetization m is 
peaked at +ms and — m^; the peak widths approach zero 
in the thermodynamic limit L — > oo, so that the average 
value mi coincides with the peak value rus (Fig. |^(a)). 

For the conserved order parameter case, the value of 
the magnetization is a constant and is same in both the 
disordered and ordered phase. One therefore needs a 
quantity that is sensitive to the difference between order 
and disorder. The simplest such quantity is the lowest 
nonzero Fourier mode of the density B 



J_|^g27riTi/L 0_ 



(3) 



where Sn denotes the average magnetization in the (d — 
l)-dimensional plane n oriented perpendicular to the x 
direction. The modulus in Eq. |^ above leads to the same 
value for all states which can be reached from each other 
by a translational shift. In the low-T ordered phase, 
Prob(Q) is expected to be a sharply peaked function, 
with peak widths vanishing in the thermodynamic limit. 
Then the mean value Qi defined by 



Qi = iimi_,oo(|Q|) (conserved) 



(4) 



serves as an order parameter. A disordered state corre- 
sponds to Qi = 0, while a perfectly ordered state with 
m = 1 in half of the system and m = — 1 in the other 
half corresponds to Qi = I/tt ~ 0.318. 
(b) Another characterization of the order is obtained 
from the asymptotic value of the 2-point spatial corre- 
lation function C{r) = (soSo+r)- At large separations r, 
C(r) is expected to decouple: 



Linir 



~,LimL^oo{SoSo+r) = {So){So+r) = ™c- (5) 



A finite value of rric indicates that the system has long- 
range order. A value rric = 1 would indicate a perfectly 
ordered pure phase without any droplets of the other 
species (like the T = state of an Ising ferromagnet), 
while rric 1 would indicate that the phase has an ad- 
mixture of droplets of the other species (like the state of 
an Ising ferromagnet for < T < Tc). 

In a finite system, C is a function only of the scaled 
variable r/L in the asymptotic scaling limit r — s- oo, L — s- 
oo (see also (d) below). An operational way to find the 
value of rric is then to read off the intercept (r/L ^ 0) in 
a plot of C versus r/L; it gives m^ in the L oo limit. 

In equilibrium systems of the type discussed above, mi 
(defined in Eq. ||) and rric coincide. 
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2. Characteristics of Fluctuations 

(c) With a conserved scalar order parameter, the low- 
T state is phase-separated, with each phase occupying 
a macroscopically large region, and separated from the 
other phase by an interface of width W . The interfacial 
region is quite distinct from either phase, and on the scale 
of system size, it is structureless and sharp. 

(d) Customarily in phase-ordered steady states, the spa- 
tial correlation function C(r) has a scaling form in Ir/Lj, 
for ^ « r « L where L is the size of the system. In 
the limit r co.L oo, \r/L\ 0, C(r) follows the 
form 0] 



C(r)«m^(l-2|r/L|) (|r/i| -> 0) 



(6) 



The origin of the linear fall in Eq. ^ is easy to understand 
in systems where phases are separated by sharp bound- 
aries on the scale of the system size, as in (c) above: a 
spatial averaging of SgSo+r produces -t-?n^ with probabil- 
ity (1 — |r/L|) (within a phase) and — with probability 
\r/L\ (across phases). The linear drop with \r/L\ implies 
that the structure factor S{k), which is the Fourier trans- 
form of C(r), is given, for large wave-vectors {kL >> 1), 
by: 



m 



{kL) 



d+l ■ 



(7) 



This form of the decay of the structure factor for scalar 
order parameters is known as the Porod law. 

It is worth remarking that the forms Eqs. ^ and |^ 
also describe the behaviour of the two-point correlation 
function in an infinite system undergoing phase order- 
ing starting from an initially disordered state. In such 
a case, L denotes the coarsening time-dependent length 
scale which is the characteristic size of an ordered do- 
main. 

(e) For usual phase-ordered systems, spatial fluctuations 
are negligible in the limit of the system size going to infin- 
ity. Hence the averages of 1-point and 2-point functions 
over an ensemble of configurations are well represented 
by a spatial average for a single configuration in a large 
system. 



B. Fluctuation-Dominated Ordering 

The phase ordering of interest in this paper occurs in 
certain types of nonequilibrium systems, and the result- 
ing steady state differs qualitatively from the ordered 
state of equilibrium systems and other types of nonequi- 
librium systems considered earlier ||^. The primary dif- 
ference lies in the effects of fluctuations. Customarily, 
fluctuations lead to large variations of the order param- 
eter which scale sublinearly with the volume, and so are 
negligible in the thermodynamic limit. Fluctuation ef- 
fects are much stronger here, and lead to variations of 



the order parameter in time, without, however, losing 
the fact of ordering. Below we discuss how the proper- 
ties (a)-(e) discussed above are modified. 

(a) Nonzero values of the averages Ml and Qi (Eqs. | 
and ^) continue to indicate the existence of order, but 
no longer provide an unequivocal measure of the order 
parameter. This is because the probability distributions 
Probijn) and Prob{Q) remain broad even in the limit 
L — > cx) (as shown schematically in Fig 0(6)). 

(b) The measure \mc\ of long-range order is nonzero, and 
its value can be found from the intercept C{\r/L\ 0). 
However, the value of rUc is, in general, quite different 
from mi. 

(c) As with usual ordered states, the regions of pure 
phases are of the order of system size L. But in contrast 
to the usual situation, there need not be a well-defined in- 
terfacial region, distinct from either phase. Rather, the 
region between the two largest phase stretches is typi- 
cally a finite fraction of the system size, and has a lot 
of structure; this region itself contains stretches of pure 
phases separated by further such regions, and the pat- 
tern repeats. Representative spin configurations {si} for 
the two cases are depicted schematically in Fig. g. This 
nested structure is consistent with a power-law distri- 
bution of cluster sizes, and thus of a critical state. The 
crucial extra feature of the FDPO state is that the largest 
clusters occupy a finite fraction of the total volume, and 
it is this which leads to a finite value of nic as in (b) 
above. Representative spin configurations |si} for the 
two cases are depicted schematically in Fig. 




FIG. 2. Schematic depiction of a linear and a cuspy decay 
of C(r) as a function of r/L, characteristic of normal phase 
ordering and FDPO respectively. Typical configurations cor- 
responding to the two cases are also shown, with 1 and — 1 
denoting the values of Si. 

(d) The ensemble-averaged spatial correlation function 
C(r) continues to show a scaling form in \r/L . However, 
in contrast to Eq. |^ it exhibits a cusp (Fig. |) at small 
values of |r/i|: 



C{r)^ml{l~b\^n. 

This implies that the scaled structure factor varies as 
S{k) 1 



(8) 



(kir 



(9) 
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with a < 1. This represents a marked deviation from the 
Porod law (Eq. |^). We will demonstrate in some cases 
that this deviation is related to the power-law distribu- 
tion of clusters in the interfacial region separating the 
domains of pure phases, as discussed in (c) above, 
(e) The spatial average of 1-point functions (m or Q) 
and the 2-point function C as a function of \r/L\ in a 
single configuration of a large system typically do not 
represent the answers obtained by averaging over an en- 
semble of configurations. This reflects the occurrence of 
macroscopic fluctuations. 



C. Fluctuating Surfaces and Sliding Particles 

Having described the general nature of fluctuation- 
dominated phase ordering, we now discuss the model sys- 
tems that we have studied and which show FDPO. We 
consider physical processes defined on a fluctuating sur- 
face with zero average slope. The surface is assumed to 
have no overhangs, and so is characterized by a single- 
valued local height variable h{x,t) at position x at time 
t as shown in Fig. ^. The evolution of the height profile 
is taken to be governed by a stochastic equation. The 
height-height correlation function has a scaling form |^] 
for large separations of space and time: 

{[Kx,t) - h{x\t')Y) ^ \x - x'^'^fi^^^). (10) 

\x — x'y 

Here / is a scaling function, and x ^-i^d z are the rough- 
ness and dynamical exponents, respectively. A common 
value of these exponents and scaling function for sev- 
eral different models of surface fluctuations indicate a 
common universality class for such models. In this pa- 
per we will study one-dimensional surfaces belonging to 
three such universality classes of surface growth. Similar 
studies of two-dimensional surfaces j^] show that similar 
fluctuation-dominated phase-ordered states arise in these 
cases as well. 



h(x,t) 



s(x,t) - - - - + + + + + + + + - - + + + + - - - - 

FIG. 3. Schematic depiction of a surface with height hix, t) 
and the coarse-grained depth function s{x,t). 

Before turning to the physical model of particles slid- 
ing on such fluctuating surfaces, we address the notion 
of phase ordering in coarse-grained depth models associ- 
ated with these surface fluctuations. In Fig. ^ we show 
the function s{x,t) which take values -1-1, —1 and de- 
pending on whether the height is below, above or at the 
same level as some reference height (h). Explicitly, we 
have s{x,t) — —sgn(h(x,t) — (h)). Different deflnitions 
of (h) define variants of the model; these are studied in 
Section II. 




Starting from initially flat surfaces, we study the coars- 
ening of of up-spin or down spin phases, which arise from 
the evolution of surface profiles. With the passage of 
time, the surface gets rougher up to some length scale 
C{t). The profile has hills and valleys, the base lengths 
of which are of the order of C(t), implying domains of 
like-valued s whose size is of the same order. Once the 
steady state is reached, there are landscape arrangements 
of the order of the system size L which occur on a time 
scale L^. However, these landscape fluctuations do not 
destroy long-range order, but cause large fluctuations in 
its value. 

Now let us turn to the problem of hard-core particles 
sliding locally downwards under gravity on these fluctu- 
ating surfaces. Figure ^ depicts the evolution of particles 
falling to the valley bottoms under gravity. When a lo- 
cal valley forms in a region (Fig. 3(a) Fig. 3(b)), 
particles in that region tend to fall in and cluster to- 
gether. The point is that particles stay together even 
when there is a small reverse fluctuation (valley — )■ hill 
as in Fig. ^(b) (c)); declustering occurs only if there 
is a rearrangement on length scales larger than the size 
of the valley. The combination of random surface fluc- 
tuations and the external force due to gravity drive the 
system towards large-scale clustering. Results of our nu- 
merical studies show that in the coarsening regime, the 
typical scale of ordering in the particle-hole system is 
comparable to the length scale over which surface rear- 
rangements take place. Further, the steady state of the 
particle system exhibits uncommonly large fluctuations, 
reflecting the existence of similar fluctuations in the un- 
derlying coarse-grained depth models of the hill-valley 
proflle. Similar effects are seen in 1-point and 2-point 
correlation functions. 



(a) 



(b) 





(c) 

FIG. 4. Depicting clustering of particles (•) in a section 
of the fluctuating surface. A surface fluctuation such as 
(a) (b) causes the particles to roll into a valley. They 
remain clustered even after a local reverse surface fluctuation 
(6) (c) occurs. 



II. FDPO IN COARSE-GRAINED DEPTH (CD) 
MODELS OF SURFACES 
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A. Surface evolution 

The dynamics of surface fluctuations can be mod- 
elled by Langevin-type equations for the height field 
h{x,t). The evolution equations for the one- 
dimensional Edwards- Wilkinson (EW) Kardar- 
Parisi-Zhang (KPZ) and Golubovic-Bruinsma-Das 
Sarma-Tamborenea (GBDT) surface fluctuations are 
respectively 



KPZ 



dh 




9t 


^ dx"^ 


dh 


d^h 


at 


^ dx'^ 


dh 








'at 


dx^ 



r]i{x,t) 



(11) 



where r]i{x,t) is a white noise with (771) — and 
{rii{x',t')r]i{x,t)) = r5{x' - x)6{t' - t), and i/i, A and K 
are constants. 

In one dimension, the EW and KPZ models can be 
simulated using lattice gas models whose large-distance 
large-time scaling properties coincide with those of the 
corresponding continuum theories. The lattice gas is 
composed of ±l-valued variables {rj_j.} on a 1 — c? lattice 
with periodic boundary conditions, where the t spins oc- 
cupy the links between sites. The values Ti-^ = +1 or 
— 1 represent the local slopes of the surface (denoted by 
/ or \, respectively). The dynamics of the interface is 
that of the single-step model , with stochastic corner 
flips involving exchange of adjacent r's; thus, /\ \/ 
with rate pi, while \/ — s- /\ with rate gi. For symmet- 
ric surface fluctuations {pi = gi), the behavior at large 
length and time scale is described by the continuum EW 
model. For pi 7^ 91, the surface evolution belongs to the 
KPZ class. Corresponding to the configuration {t^^i} 

we have the height profile {hi} with hi — J2i<j<i'''j-^- 

For simulating a surface fluctuating via a GBDT pro- 
cess, we used a solid-on-solid model with depositing par- 
ticles piling up on top of each other. The height hi at site 
i is the height of the pile of particles at that site. Dur- 
ing each micro-step a particle is deposited randomly on a 
site i. If the new height hi at i, is greater than hi-i and 
/li+i, then with equal probability (= 1/3) three things 
are attempted — the deposited particle can remain at 
site i, or can move to the neighboring sites i — 1 or i + 1. 
It actually completes the left or right move only if there 
is an increase in the coordination-ordination number of 
the particles |6|. 



B. Definitions of the CD models 

Let us imagine a process of coarse-graining which elim- 
inates fine fluctuations of the height profile, and replaces 
the height field hi at site i by a variable Si which is -1-1, 
-1 or depending on whether the surface profile at site i 



is below, above or exactly coincident with a certain refer- 
ence level, which is the same at all i. The aim is to have 
a coarse-grained construction of locations of large valleys 
and hills. Our procedure depends on the choice of the 
reference level, and we have explored three choices (the 
GDI, CD2 and CDS models) which are discussed below. 

In model GDI, the reference level is set by the initial 
condition, which corresponds to an initially flat inter- 
face: h{x,t = 0) = 0. The coarse-grained depth function 
is then 

s{x,t) ^ -sgn[h{x,t)]. (12) 

With the passage of time, the surface becomes rougher, 
so that h{x, t) develops hills and valleys with respect to 
the level. As the base lengths of the hills and valleys 
grow in size, there is a growth of the domains of the vari- 
able s{x,t). We are able to characterize the coarsening 
behaviour of this model analytically in some cases. 

In a finite system, at long enough times the surface 
moves arbitrarily far away from its initial location. Thus 
the steady state of the GDI model is trivial — all Si are 1, 
or all are -1, with probability one. This clearly happens 
because the reference level in the GDI model is fixed in 
space. This leads us to examine models (GD2 and GD3) 
where the reference level moves along with the surface, 
so that we may expect nontrivial steady state properties. 

In model GD2, the coarse-grained depth function 

Si = -sgn[hi] (13) 

where hi ~ Si<j"<i''j-i defined in the earlier sec- 
tion. Note that at all times t, the origin is pinned so 
that /ii=o = 0. The height function of the continuum 
version of the GD2 model is related to that of GDI as: 
hf"^{t) = hf"\t)-h^"\t). The function is +1, -1 
or accordingly as the height hi at site i is below, above 
or at the zero level. A stretch of like s^'s = -1-1 represents 
a valley with respect to the zero level. The time evolution 
of the GD2 model variables {si} is induced by the under- 
lying dynamics of the bond variables {Tj_i} defined in 
the previous subsection. This model was studied by us 

Finally model GD3, is defined as follows: hi is con- 
structed from r's exactly as described for the GD2 model, 
but then one defines 



-sgn[hi - {h{t))] 



(14) 



where {h{t)) — (1/L) hi{t) is the instantaneous av- 

erage height which fluctuates with time. This definition 
was used earlier by Kim et al jl^ ] who were studying 
domain growth in an evolving KPZ surface. 

Each of the GD models defined above has its own mer- 
its and limitations. We will see below that the GDI 
model proves to be analytically tractable (for Gaussian 
surface fluctuations) in the coarsening regime, while for 
the GD2 model several exact results can be derived in 
the steady state. Of the three models, the GD3 model 
most resembles the model of sliding hard-core particles 
on the surface that is studied in Section |^ below. 
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C. Coarsening in the CD models 

1. Analytical results for the CDl model 

In this section, our primary focus is on coarsening 
properties of a class of CDl models. To this end, we 
will focus on the equal time correlation function 

C(x,t) = {a{Q,t)a{x,t)) = {sgn[h[0,t)]sgn[h{x,t)]) . 

(15) 

We consider only linear interfaces evolving from a flat 
initial condition h{x, 0) = according to the Langevin 
equation. 



dh 
'dt 



(16) 



where r/(x, t) is a Gaussian white noise with {ri{x, <)) — 
and {r]{x,t)r]{x' ,t')) — 6{x — x')5{t — t'). The dyanmic 
exponent z specifies the relaxation mechanism. For ex- 
ample, z — 2 corresponds to Edwards- Wilkinson (EW) 
interface and z — A corresponds to Golubovic-Bruinsma- 
Das Sarma-Tamborenea (GBDT) interface. Since tj(x, t) 
is a Gaussian noise and the evolution equation ( p^ ) is 
linear, the height field h{x,t) is a Gaussian process. For 
Gaussian processes, it is straightforward to evaluate the 
correlation function in Eq. (Qq) exactly and one finds. 



C{x,t) = -sin^i [H{x,t)] 

TT 

where H(x,t) is given by, 

(h{0,t)h{x,t)) 



H{x,t) 



(17) 



(18) 



Now the normalized height correlation function H{x,t) 
can be easily computed for linear interfaces evolving 
via Eq. ( [l6|) by taking the Fourier transform of Eq. 
(p^). From Eq. (^6|), assuming flat initial condition, the 
Fourier transform, {h{k,t)h{—k,t)) is given exactly by. 



{h{k,t)h{-k,t)) 



(1 



-2|fc|^t 



2\k\' 



Inverting this Fourier transform we get 

(z-1 



H{x,t) 



2i-ir(i 



-F 



t- 



where the scaling function F{y) is given by, 



1 



■ cos{yu)du. 



(19) 



(20) 



(21) 



Using this exact expression of H{x,t) in Eq. (Il7|), we 
get the exact correlation function for arbitrary linear in- 
terface model parametrized by the dynamic exponent z. 



It is also evident that C{x, t) is a single function of the 
scaled distance, y = xt 

The small distance behavior of the scaling function can 
be easily derived from the small argument asymptotics of 
the integral in Eq. (|2l|). Let us first consider the EW in- 
terface with z = 2. In this case the integral in Eq. ([2l| ) 
can be done (by putting a factor w in the exponential. 



i.e., writing e 



and then differentiating with respect 



to w and then integrating back with respect to w upto 
w = \) and we get. 



(22) 



A change of variable, x'^/Swt — y gives a more compact 
expression, 



H{x,t) 



4V2t J^2/st 



e-yy~^^^dy. 



(23) 



Integration by parts yields the desired short distance be- 
haviour. 



H{x,t)^l- J-\x\ 



(24) 



Putting this back in Eq. (10) and expanding the arc- 
sine, we get. 



C{x,t) = 1 



xt' 



+ . . 



(25) 



Thus the correlation function has a square-root cusp at 
the origin for the z — 2 CDl model. One can similarly 
do the small distance analysis for arbitrary z > 1. We 
find that for general z, 



C{x,t) = l-a|a;t-i/"|" + 



(26) 



where a is a z-dependent constant and the cusp exponent 
a is given by. 



a — {z 
a 1 



-l)/2 for 
for z > 3. 



z < 3 



(27) 



For z — 3, we find additional logarithmic corrections. 



C{x,t) = l-a\y\^ 



y\ 



(28) 



where y — xt~^/'^ . 

Thus our exact results indicate that z = Zc = 3 is a 
critical value. For z > 3, one recovers the linear cusp 
in the correlation function at short distances (and hence 
Porod's law) indicating sharp interfaces between domains 
as in the usual phase ordering systems. But for z < 3, 
one gets a z-dependent cusp exponent signalling anoma- 
lous phase ordering dominated by strong fluctuations and 
a significant deviation from Porod's law. The value Zc is 
the one across which a morphological transition has been 
shown to occur in Gaussian surfaces p9[ |, in the context 
of spatial persistence of fluctuating surfaces. 
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2. Numerical results for the CDS model 



Unlike the GDI model, we have not been able to an- 
alytically characterize the coarsening properties of the 
CD2 or CDS models, in which the reference level moves 
with time. However the coarsening properties in both 
CD2 and CDS models can be studied numerically. Be- 
low we present the numerical results for the equal time 
correlation function C for the CDS model in three dif- 
ferent cases where the underlying surface is evolving re- 
spectively by the EW, KPZ and GBDT dynamics. The 
initial condition chosen was i = 1 at odd bonds and 

J 2 

— 1 at even bond locations, ensuring that the height pro- 
file was globally flat. We used a lattice with a number 
L = 409600 of bonds and equal number of sites. At time 
t > correlations gradually develop as the s-spin do- 
mains grow. In Figs. ||, ^, and ^ we show the data for 
C as a function of t (insets of the respective figures), and 
how they collapse on to a single curve Cs in each case, on 
scaling r by a ^-dependent length scale C{t). For each of 
the three cases, we see that C{t) ~ t^^^, where the dy- 
namical exponent z = 2, S/2, and 4, respectively for the 
EW, KPZ and GBDT surfaces. Notice that the scahng 
curves for EW and KPZ surfaces have a cusp at small 
values of the argument r/C, and the cusp exponent (Eq. 
|, 1) a ~ 0.5 for both. For the GBDT surface there is 
no cusp, and a ~ 1.0. We note that these results for the 
CDS model are consistent with the analytical results in 
Eq. @ of the GDI model. 

The fact that the correlation function has a scaling 
form in r/C{t), with a nonzero intercept implies that at 
infinite time the system would reach an ordered steady 
state, as the value of C at any fixed r (no matter how 
large) approaches the value of the intercept at large 
enough time. The intercepts of all the three curves in 
Figs. ^, ^, and ^ have the value 1 implying that nic = 1 
for the CDS model. 
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FIG. 5. The data shown in the inset for C(r, t) for the CDS 
model of the EW surface at different times t = 400x2" (with 
n = 0,...,6), is seen to collapse when r is scaled by £,{t) ~ t^^^. 
The cusp in the scahng function at small argument is charac- 
terized by a ~ 0.5. 
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FIG. 6. The data shown in the inset for C(r, t) for the CDS 
model of the KPZ surface at different times t = 400x2" (with 
n = 0,...,6), is seen to collapse when r is scaled by £(t) ^ t}^^ . 
The cusp in the scaling function at small argument is charac- 
terized by a ~ 0.5. 
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FIG. 7. The data shown in the inset for C(r, t) for the CDS 
model of the GBDT surface at different times t = 400x2" 
(with n = 0,...,6), is seen to collapse when r is scaled by 
£(t) ~ t^''*. The behaviour of the scaling function at small 
argument is characterized by a ~ 1.0. 



D. Steady state of the CD models 

In a finite system, as time passes the surface diffuses 
away from its i = location. As discussed above, this 
leads to a trivial steady state in the GDI models, corre- 
sponding to all Si — 1 (or all s; — —1) with probability 
one. We need the reference level to keep up with the sur- 
face in order to probe the steady state aspects of coarse- 
grained surface fluctuations. This is accomplished in the 
CD2 and CDS models. 

In both the CD2 and CDS models we will see below 
that the cluster size distribution of the Si variables varies 
as a power law ~ in the steady state. The order pa- 
rameters have a broad distribution and the scaled 2-point 
function has a cusp. 

It is well known that for both EW and KPZ surfaces 
in 1 — d, the steady states have random local slopes [|| , 
i.e. the steady state probability distribution of the height 
profile is 



P{{K\) = P„e 



(29) 
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This leads to a mapping of each surface configuration in 
the CD2 and CDS models to a random walk (RW) tra- 
jectory. The correspondence is as follows: Ti_i/2 = +1 or 
— 1 can be interpreted as the rightward or leftward RW 
step at the i'th time instant. Then in the CD2 model, 
Si = 1,-1 or depending on whether the walker is to 
the right, to the left or at the origin after the z'th step. 
In the CDS model, the reference point for demarcating 
left (si = 1) and right {si = —1) is the average of dis- 
placements (heights), and can be fixed only after the full 
trajectory is specified; then with respect to (/i), the value 
of the position of the walker at every I'th instant gets 
specified and hence also the Si spins. 



1. Power Law Distribution of Cluster Sizes 

For the CD2 model with EW or KPZ dynamics, ex- 
act results for different properties in the steady state 
can be derived, because the surface profiles map on to 
random walks. Periodic boundary conditions imply that 
the RW starts at time from the origin and comes back 
to the origin after L time steps. Evidently, the lengths 
of clusters of s = 1 spins (or s — ~l spins) represent 
times between successive returns to the origin. Thus 
P{1), the probability distribution of the cluster sizes 
for the CD2 model is exactly the well-known distribution 
(« l/vzTiT^e^^/^^')) for RW return times to the origin, 
which behaves as ^ l"^^'^ (for large /) with a cutoff at 
I = L. Thus 6* = S/2 in this model. 

For the CDS model, the variable reference point makes 
it difficult to make exact statements, but we expect that 
the cluster size distribution at large lengths I will still 
be given as The numerically determined P{iys for 

the CD2 and CDS models are plotted in Fig. ||, and they 
show the expected power-law decay. 

We note that the power-law distribution of the inter- 
vals between successive returns is related to the spa- 
tial persistence of fluctuating interfaces For lin- 
ear interfaces evolving via Eq. (|l^), the corresponding 
spatial persistence exponents were computed recently in 
Ref. ||l9[. From these results we infer that for CD mod- 
els with GBDT dynamics, the cluster size distribution 
in the steady state also has a power law distribution, 
P{1) - l-^/'^ in (1 + l)-dimensions. 



FIG. 8. P{1) against I for up-spin clusters decays as I 
in both the CDS (•) and CD2 (empty triangles) models. We 
used L = 2048. 



2. Order Parameter Distribution 



The distributions of the order parameters for each of 
the CD2 and CDS models are broad. For the CD2 model, 
an appropriate (nonconserved) order parameter is aver- 
age value mi of modulus of m = x S (^^^ ^1- H) ' 
which for the RW represents the excess time a walker 
spends on one side of the origin over the other side. In 
order to respect periodic boundary conditions, we need to 
restrict the ensemble of RWs to those which return to the 
origin after L steps. The full probability distribution of 
m over this ensemble is known from the equidistribution 
theorem on sojourn times of a RW fl^ : 



Prob{m) = 1/2 



m e [-1, 1] 



(SO) 



i.e. every allowed value of m is equally likely. This im- 
plies (|m|) = 1/2 and ((m^) - {\m\fy/^ = l/\/T2. 

For the CDS model, most often half of the surface pro- 
file is above the average height level and half below it. 
As a consequence, we find numerically that the distri- 
bution of cluster sizes P{1) decays sharply beyond L/2. 
This resembles the sliding hard-core particles and hence 
the conserved order parameter Qi is more suitable to de- 
scribe the ordering in this model than mi. We monitor 
the average value Qi of Q* = j^'^'^''^^ Pj\ (Eq. |), 
where pj = {l+Sj)/2. This order parameter has a value 
for a disordered configuration and a value I/tt « 0.S18 for 
a fully phase separated configuration with two domains 
of -I- and — spins, each of length L/2. The numerical 
value of the distribution P{Q*) of Q* is shown in Fig. 
^, and the average value Qi in the limit of large system 
size numerically approaches the value 0.22. It is apparent 
from Fig. ^ that P{Q*) is broad, and is larger for larger 
Q* . The width, which remains finite in the thermody- 
namic limit, signifies that large-scale fluctuations occur 
frequently in the system. 



0.08 



0.01 



0.0001 



10 100 1000 

1 



0.06 ^ 



O 0.04 



0.1 



0.2 



0.3 



Q* 



FIG. 9. Probability distribution P{Q*) in the steady state 
of the CDS model The mean value is Qi ~ 0.22. 
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3. Correlation Functions 



A. Correlation Functions through the Independent 
Interval Approximation 



Finally we turn to the 2-point spatial correlation func- 
tions in the steady state of CD2 and CDS models. The 
growing length scale C{t ) as t — > oo is limited by the sys- 
tem size L. In Fig. we show the scaling of data for 
C{r) in the steady state as a function of r/L for an EW 
surface. Both the curves show a cusp at small values of 
r/L, with cusp exponent a ~ 0.5. 



We now show analytically within the Independent In- 
terval Approximation (IIA) that the cusp exponent 
a and the power law exponent 9 are related. Within this 
scheme, the joint probability of having n successive inter- 
vals is treated as the product of the distribution of single 
intervals. In our case, the intervals are successive clusters 
of particles and holes, which occur with probability P{1). 
Defining the Laplace transform P{s) — dle^''^ P[l), 
and C{s) analogously, we have Q 
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FIG. 10. Steady state C(r) collapses onto a single curves 
when plotted against r/L for both the CD2 and CDS mod- 
els.. The scaling function shows a cusp at small values of the 
argument, with a ~ 0.5 for both models. We used L — 64, 
128, 256 and 512. 



Since successive RW returns to the origin are indepen- 
dent events, the calculation in section III A below based 



on independence of intervals, is exact for the CD2 model. 
Thus Eq. |3^ holds, and we conclude that the correlation 
function cusp exponent a = 1/2 exactly for the CD2 
model. This also implies the result that a = 1/2 even in 
the coarsening regime for the CD2 model with EW and 
KPZ surfaces. This is because at any time t, regions of 
a coarsening system which get equilibrated are of length 
~ C{t) << L. Now the correlation function C{r,t) is ob- 
tained by spatial averaging over the system, and hence 
equivalently averaging over an ensemble of several steady 
state configurations of subsystem size ^ C{t). Thus the 
exact result for a in the steady state carries over to the 
coarsening regime. 



s(l ~ sC{s)) 



2 l-P(s) 
(0 1 + P(s) 



(31) 



where (/) is the mean cluster size. In usual applica- 
tions of the IIA, the interval distribution P{1) has a 
finite first moment (/) independent of L. But that is 
not the case here, as P{1) decays as a slow power law 
P{1) - l-^eiL-l) for / >> 1, with the function Q{L~l) 
denoting that the largest possible value of Z is i. This 
implies that (/) « aL^~^ for large enough L. Con- 
sidering s in the range 1/L << s << 1, we may ex- 
pand P{s) « 1 — bs^~^; then to leading order, the right 



hand side of Eq. (31) becomes bs^ ^/aL^ ^, implying 



C(s) « 1/s - b/{aL^-\s^-''). This leads to 



Cir) 



\2-0 



aT{3-0)'L' 



(32) 



This has the same scaling form as Eq. (||) . Matching the 
cusp singularity in Eqs. (||) and (|3|), we get 



a = 2 



(IIA). 



(33) 



We recall (see Section 11 D ) that the assumption of inde- 
pendent intervals which underlies the IIA in fact holds 
exactly for the CD2 model, and Eq. |3^ implies that 
a = 1/2 in the steady state and the coarsening regime for 
the CD2 model. For other models like the CD3 model, 
or the sliding particle models we will encounter in the 
subsequent sections, the IIA gives insight into the origin 
of the cusp from the power laws, although it is not exact. 



III. UNDERSTANDING FDPO IN CD MODELS 



We have seen in the previous section that the distribu- 
tion of like-spin clusters follows a slow power law decay 
in the CD models. We will demonstrate below that on 
the basis of this power law, we may understand both (i) 
the cusp in the 2-point function, as well as (ii) ordered 
phases which occupy a finite fraction of system size. 



B. Extremal Clusters and Ordered Phases 

We now turn to our claim (ii), that the very same dis- 
tribution which gives rise to power-law distributed broad 
boundaries with a collection of small clusters, also gives 
rise to large clusters of size '--^ L of 'up' or 'down' spins, 
which form the pure phases. For the CD2 model, we 
numerically studied the sizes of the largest cluster li for 
systems of different sizes L; we show them in Fig. |l|. 
The full distribution P{li) scales as a function of h/L. 
The average value is (^i) ~ 0.48L. We also find a similar 
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scaling of the distribution Pih), for the second largest 
clusters of size I2, and {I2) — 0.16L (see Fig. |ll|) . 



Some understanding of the fact that the size of largest 
clusters are of order L can be reached by considering the 
statistics of extreme values. Applied to our case, if N 
cluster lengths are drawn at random from a distribution 
of lengths given by P{1) ^ {0 — 1)//^, then the probabil- 
ity distribution L m [x) that the largest cluster is of length 
X goes as ~ N x~^ exp{— N x~'^^~^^) pO]. The latter dis- 
tribution peaks at a; = Xmax ]Sf^'^~^\ In the CD2 
problem 6 = 3/2. Now, in a system of length L we have 
on an average a/L clusters. If we make the approximate 
replacement of N by this average number VL, we imme- 
diately get Xmax ~ L. This explains how, although the 
average cluster sizes are of order L^/^, there are always 
clusters with sizes of order L. This is reminiscent of the 
behaviour of the largest loops in a random walk . 




FIG. 12. Estimates of the magnetization mi (squares), 
m2 (triangles) and rh^ (•) from the largest few clusters plot- 
ted against the total magnetization for m, for 1000 different 
configurations. The convergence towards the line of slope 
unity shows that a few large clusters account for the major 
contribution to m. 



IV. HARD-CORE PARTICLES SLIDING ON 
FLUCTUATING SURFACES 



A. The Sliding Particle (SP) Model 



'i.. 
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1/L 

FIG. 11. Probability distributions P(/i) (the curve on the 
right) and Pih) of the largest (length h) and second largest 
(length h) clusters in the steady state of CD2 model are seen 
to collapse when plotted against scaled lengths h/L and h/L, 
respectively. The sizes used are L = 512, 1024, 2048 and 4096. 



Further we found the contribution to magnetization 
coming from the largest clusters in the system and com- 
pared them with the total magnetization of the system, 
configuration by configuration. In Fig. |l^, we show 
scatter plots of mi which is the magnetization obtained 
from summing the spins of the largest cluster, 7712 which 
is obtained by summing spins of largest and the sec- 
ond largest cluster, and 7713 by summing those down to 
the third largest cluster against the total magnetization 
m — {i/L)J2^i- The convergence of the scatter plots 
towards the 45° line, shows that the few largest clusters 
give a major contribution to the magnetization of the 
system. Each of these large clusters is a pure phase with 
magnetization 1, and thus gives rise to rric = 1 in the 
curves in Fig nO. 



In this section we consider the physical problem of 
hard-core particles sliding locally downwards on the fluc- 
tuating surfaces discussed in the previous sections. We 
find that the downward gravitational force combined with 
local surface fluctuations lead to large scale clustering of 
the hard-core particles. The phase-separated state which 
arises mirrors the hill- valley profile of the underlying sur- 
face. For example, the particles on EW and KPZ surfaces 
show FDPO with the cluster distribution, 1-point fimc- 
tion, and 2-point function behaving as in their CD model 
counterparts. On the other hand, particles on the GBDT 
surface show conventional ordering. 

Let us first define a sliding particle (SP) model on a 
one-dimensional lattice. This is a lattice model whose 
behaviour resembles that depicted in Fig. ^. The parti- 
cles are represented by ±l-valued Ising variables {cri} on 
a one-dimensional lattice with periodic boundary condi- 
tions, where a spins occupy lattice sites. The Tj_i vari- 
ables occupy the bond locations and represent the sur- 
face degrees of freedom as described in section II B for 
the CD2 model, and their dynamics involves indepen- 
dent evolution via rates pi and qi as discussed earlier. 
For the particles, pi = ^(1 -I- cr^) represents the occupa- 
tion of site i. A particle and a hole on adjacent sites 
-\- 1) exchange with rates that depend on the inter- 
vening local slope t^-I] thus the moves •\o o\» and 
o/» —>■ •/o occur at ratep2, while the inverse moves occur 
with rate (72 < P2- The asymmetry of the rates reflects 
the fact that it is easier to move downwards along the 
gravitational field. For most of our studies we consider 
the strong-field (92 = 0) hmit for the particle system. We 
set p2 = Pi- The dynamics conserves a and X] we 
work in the sector where both vanish. This corresponds 
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to a 1/2 filled system of particles on a surface with zero 
average tilt. For the EW surface, we took pi = qi, while 
for the KPZ surface we took pi = 1 and qi = 0. In Sec- 
tion 0, we discuss departures from these conditions and 
explore the robustness of FDPO to these changes. 

For the GBDT sur face , the evolution of which was 
described in Section [I A , a chosen particle moves to 
its right or left with equal probability (= 1/2) if there 
is locally a non-increasing height gradient. Thus again 
^2 = 0. The rate of update of the particles is same as 
that of the surface. 

The problem can be specified at a coarse-grained meso- 
scopic level by the continuum equations for the density 
field p{x, t) corresponding to the discrete variable pi for 
the particles. Since the particle density is conserved, 
the starting point is the continuity equation dp/dt = 
—dJ{x,t)/dx, where J is the local current. Under the 
hydrodynamic assumption, the systematic part of the 
above current is —pdh/dx, since for viscous dynamics, 
the speed is proportional to the local field, in this case 
the local gradient of height. Moreover there is a diffusive 
part —V2dpldx which is driven by local density inhomo- 
geneities, and a noisy part 772(2;, which arises from the 
stochasticity. The noise 772 is a Gaussian white noise. 
The total density can be written a,s p = po + p, where po 
is the average density and p is the fluctuating part. This 
implies finally that the density fluctuation p evolves via 
the following equation: 



dp 

dt 



^ dx dx' 
dri2{x,t) 



+ (1 - 2p, - m±.)i-) + (1 - 2po)p(Tr2 ) 



' dx^ 



dx 



(34) 



Using the well-known mapping in 1 — d between the den- 
sity field p and the height field h of the corresponding 
interface problem jl^, one has the relation p ~ dh/dx. 
This implies from Eq. ^ that the lowest order term 
in the evolution equation of h is proportional to dh/dx. 
This linear first-order gradient term is the result of the 
gravitational field which acts on the particles. The evo- 
lution of the field h{x, t) is given by Eq. Thus a con- 
tinuum approach to the problem of the sliding particles 
requires analysis of the semi-autonomous set of nonlinear 
equations ^ and ^ as one of the fields evolves indepen- 
dently but influences the evolution of the other. The 
problem belongs to the general class of semiautonomous 
systems, like the advection of a passive scalar in a fluid 
system JTsf . 

The SP model is a special case of the Lahiri- 
Ramaswamy (LR) model [ p"7|Jl8| of driven lattices such 
as sedimenting colloidal crystals. The general LR model 
has two-way linear couplings between the p and h fields, 
and its phase diagram has recently been discussed in [ p2[ . 
The SP model of interest here has autonomous evolu- 
tion of the {h{x)}, and corresponds to the LR critical 



line which separates a wave- carrying phase p3[ from a 
strongly phase separated state ||l^. Further, in a model 
of growing binary films considered in [ p4[ , in the limit 
where the height profile evolves independently, the prob- 
lem gets mapped to noninteracting domain walls (if an- 
nihilation is neglected) rolling down slopes of indepen- 
dently growing surfaces. The latter problem becomes 
similar to ours, on thinking of the domain walls as parti- 
cles. But the fact that they are noninteracting in contrast 
to the hard-core particles may introduce other physical 
effects into the problem. 



B. Coarsening in SP model 



We start with a surface in steady state, and allow an 
initially randomly arranged assembly of sliding particles 
to evolve on it. In an initial short-time relaxation, par- 
ticles slide down to the bottom of local minima. After 
this, the density distribution evolves owing to the rear- 
rangement of the stochastically evolving surface, whose 
local slopes guide particle motion. We found in nu- 
merical simulations that the surface fluctuations actu- 
ally drive the system towards large scale clustering of 
particles. This can be seen as follows. After time t, 
the base lengths of coarse-grained valleys of length t^^^ 
would have overturned, where z is the dynamical expo- 
nent of the surface. We thus expect that the latter length 
scale sets the scale of particle clustering at time t. To 
test this we monitored the equal time correlation func- 
tion C{r,t) = {ao{t)(To+r{t)) by Monte-Carlo simulation. 
We found that it has a scaling form 



C = f{r/C[t)) 



with C ^ t 



l/z 



(35) 



in accord with the arguments given above. The data for 
C(r, t) for the particles on EW, KPZ and GBDT surfaces 
are shown to collapse in Figs. 1^ and |l5|, respectively. 
Evidently, Eq. |3| holds quite well for all three surfaces, 
despite the widely different values of z for the three. The 
onset of scaling will be discussed further in section 0, 
where we discuss the effect of varying the ratio of rates 
of relative updates of the particles and the surface. 
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FIG. 13. The data shown in the inset for C{r,t) for the SP 
model with an EW surface at different times t — 400 x 2" (with 
n = 0,...,6), is seen to collapse when scaled by £{t) ~ t^^^. 
For small arguments, the scaling function has a cusp with 
a ~ 0.5. 



in the latter case, corresponding to a real space decay 
with a cusp exponent a. A nonzero value of Co implies 
that nic 7^ 1, as nic is given by -\/l — Co- This indicates 
that the particle-rich phase has some holes and vice versa. 




r/t2/3 

FIG. 14. The data shown in the inset for C{r,t) for the SP 
model with a KPZ surface at different times t — 400x2" (with 
n = 0,...,6), is seen to collapse when scaled by C{t) ~ t^^^ . 
For small arguments, the scaling function has a cusp with 
a ~ 0.25. 




r/ti/4 



FIG. 15. The data shown in the inset for C{r,t) for the SP 
model with a GBDT surface at different times t = 400x2" 
(with n — 0,...,6), is seen to collapse when scaled by 
£{t) ~ t^''^. The scaling function has no cusp and a ~ 1.0. 

To determine the short distance behavior of the decay 
of C as a function of r/£{t), we evaluated the structure 
factor S{k) for C. For any finite C{t), we may write 



C^Coir)+Cs{r/C), 



(36) 



where Co{r) is the analytic part which decays over small 
distances r, while Cs is the nonanalytic part which scales 
as a function of r/C. We are primarily interested in Cs, 
and so need to subtract the appropriate Co from C. In 
terms of the scaled variable y — r / C, Co contributes only 
to y = 0, in the limit £ — > oo. In that limit we write 
C{y) = Csiy) -\-Co6y.i), and determine Co by seeing which 
value gives the longest power-law stretch for S/C, as 
judged by eye. In Fig. |l^ we show S for a late time, 
obtained without any subtraction and after subtraction 
of CoSyfl with Cs(0) = 0.71. The power law decay as 
~ l/(fc£)"+^, stretches over a substantially larger range 
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FIG. 16. The scaled structure factor S{k)/C versus k£, for 
the SP model with EW surface, with (•) and without (o) 
subtraction of the analytic part. 
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FIG. 17. The scaled structure factor S/C is plotted against 
kC, corresponding to the curves for n = 6 in the insets of 
Figs. ^ |l4| and |l|. The slopes at large kC for KPZ, EW 
and GBDT are —1.25, —1.5 and —2, respectively. 

In Fig. ^ we show, S corresponding to the three dif- 
ferent surfaces at t — 400 x 2^. We find that for the EW 
surface a ~ 0.5, for the KPZ surface it is ~ 0.25, and for 
the GBDT surface it is ~ 1.0. Thus there is a deviation 
from the Porod law behavior for the EW and KPZ sur- 
face fluctuations, and no such deviation for the GBDT 
surface. In all three cases, we see that the behaviour of 
the 2-point functions in the particle system resembles the 
corresponding correlation functions of the CD model for 
the underlying surface. In the KPZ case, the value of the 
exponent a ~ 0.25 is different from its value a = 1/2 
in the CD model counterpart. For the EW and GBDT 
surfaces, the values of a are ~ 0.5 and 1.0 respectively 
as in the corresponding CD models. 

Since the SP model corresponding to the GBDT sur- 
face does not exhibit anomalous behaviour of the scaled 
two-point correlation function which is a signature of 
FDPO, we do not consider it further in our subsequent 
discussion of the steady state. 
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C. steady state of the SP model 

We first study 1-point functions in order to charac- 
terize the steady state. As the system phase separates, a 
suitable quantity to study is the magnitude of the Fourier 
components of the density profile 

Qik) = \^J2^"''n,\, k = ^. (37) 

where rij = (1 + (Tj)/2 and m = 1, L — 1. A signature 
of an ordered state is that in the thermodynamic limit, 
the average values {Q{k)) go to zero for all k, except at 
k ^ 0. We monitored these averages for the system of 
sliding particles, with the average (• • •) performed over 
the ensemble of steady state configurations. In Figs, 
and |l^ we show the values of {Q{k)) as a function of k 
for various system sizes L, for the EW and KPZ surfaces 
respectively. In both the cases, for all fc ^ the value of 
{Q{k)) falls with increasing L indicating that {Q{k)) 
in the thermodynamic limit, for any fixed, finite k. But 
for k = 2tt/L, we see that the value of {Q{k = ^)) ap- 
proaches a constant. The sharpening of the curves near 
/c — > implies an ordered steady state. 



0.16 
0.12 

O* 0.08 "J . • . 
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k 

FIG. 18. {Q{k)) plotted as function of fc = 27rm/L, for 
different system sizes L = 32, 64, 128 and 256, for an EW 
surface. 



to TO = 1) as a measure of the extent of phase separa- 
tion. We have used Qi = {Q*) as the order parameter 
earlier also for the CD3 model, and note that it has been 
used earlier in other studies of phase-separated systems 
[Q. Here we find that Qi ~ 0.18 and 0.16 for particles 
on the EW and and KPZ surfaces respectively. The lat- 
ter values being nonzero indicates that the steady state 
is ordered. At the same time, the values being less than 
0.318 indicates that the states deviates substantially from 
a phase separated state with two completely ordered do- 
mains. To have a full characterization of the fluctua- 
tions which dominate the ordered state, one should actu- 
ally evaluate the probability distributions of all the Q's, 

e.g. Q* = Q{^), g(2) = Q(3) = 

We show below (in Fig. gy) one of these distributions, 
namely that of Q* = (5(^), for an EW surface. We find 
that the distribution P{Q*) remain broad (with root- 
mean-square deviation being ~ 0.07) even as L ^ oo, in- 
dicating again the dominance of large scale fluctuations. 



0.05 I . . . — — . . . 1 

0.04 - 

0.03 • _ - 

P(Q*) 

0.02 - . - 

0.01 ■ • • ' - 

t • • • • ' • 

0.1 0.2 0.3 

Q* 

FIG. 20. Numerically determined probability distribution 
P{Q*) of the order parameter Q* , obtained for the SP model 
with EW surface fluctuations, in steady state. 
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FIG. 19. {Q{k)) plotted as function of A; = 27rm/L for dif- 
ferent system sizes L = 32, 64, 128 and 256, for a KPZ surface. 



The above behavior of {Q{k)) as a function of k sug- 
gests that we take the value Q* = Q(^), (corresponding 



It is instructive to monitor the variation of Q* as a 
function of time t, for different system sizes. For an EW 
surface (Fig. ^) the value of Q* shows strong excursions 
about its average value, consistent with the broad distri- 
bution shown in Fig. |2^. The temporal separation pe- 
riod of these fluctuations of the order parameter increases 
roughly as ~ L^, but their amplitude is independent of L. 
Consequently P{Q*) approaches an L-independent form 
a& L ^ oo. A temporally oscillatory order parameter also 
has been found earlier in a model for comparative learn- 
ing [ p5[ . However the temporal behaviour in our case is 
quite different from the almost periodic fluctuation in the 
latter model, as the Fourier spectrum of the time series 
in Q*{t) in our case follows a broad power-law. We have 
not pursued a detailed study of the temporal behaviour 
any further. 
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L=128 

Q* 
L=64 




t 10000 20000 

FIG. 21. Variation of Q* with time t, for different system 
sizes L = 32, 64, 128 and 256, showing that the separation be- 
tween the fluctuations of the order parameter increases with 
L, but that the ampUtude does not vary much. 




FIG. 22. Variation of Q* = Q(l) (solid thick line), Q{2) 
(solid thin line) and Q(3) (broken line) are shown as a func- 
tion of time to show that a decrease in value of one is ac- 
companied by an increase in the others, indicating that one 
large cluster may break up into a few large ones, in the steady 
state. The system size was L — 128. 

The fluctuation of Q* in Fig. ^ gives rise to an in- 
teresting question: Does the system become disordered 
and lose the phase ordering property when the value Q* 
falls to low values? The answer is no, as is very clearly 
brought out in Fig. H in which Q* = Q(l), Q{2), Q(3) 
have been plotted simultaneously as a function of time t 
for a single evolution of the system. We observe that a 
dip in Q* is accompanied by a simultaneous rise in the 
value of either Q{2) or Q(3). This implies that whenever 
the system loses a single large cluster (making Q* small) 
either two or three such clusters appear in its place (mak- 
ing the values of Q{2) and Q{3) go up). Thus the system 
remains far from the disordered state, and always has a 
few large particle clusters which are of macroscopic size 
~ L. A numerical study showed that the average size of 
the largest particle cluster ~ 0.14i. 

We have seen above that in the SP models, the order 
parameter has a broad distribution just as in their CD 
model counterparts. We observe further that the parti- 
cle and hole cluster size distributions in the steady state 
of the SP model decay as a power-law: P{1) ^ . In 
Fig. for the EW surface, we find that the particle (de- 
noted by symbols) and hole (denoted by lines) distribu- 



tions coincide, with 9 ~ 1.8. By contrast. Fig. |2J for the 
KPZ surface shows that the particle and hole distribu- 
tions are not identical. This is because with asymmetric 
rates (pi ^ qi), the surface has an overall motion in one 
direction, such that the downward motion of the parti- 
cles and the upward motion of the holes, due to gravity, 
are no longer symmetrical. We checked that the distribu- 
tions for particles and holes get interchanged if the rates 
Pi and qi are interchanged. The exponent for the decay 
of both the particle and hole distributions is ~ 1.85. 




FIG. 23. P{1) vs. I for clusters of particles (symbols) and 
holes (lines) in the SP model with an EW surface, for differ- 
ent system sizes L — 256, 512, 1024 and 2048. P{1) decays as 
a power law with ~ 1.8. The inset shows collapsed data of 
steady state C(r) for L — 64, 128, 256 and 512 as a function 
of r/L; the scaling function has a cusp with a ~ 0.5. 
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1 

FIG. 24. P{1) Vs. I for clusters of particles (symbols) and 
holes (lines) in the SP model tha KPZ surface, for different 
system sizes L = 256, 512, 1024 and 2048. The data shows 
the existence of a particle-hole asymmetry. A power-law with 
9 ~ 1.85 has been shown along with the curves as a guide to 
the eye. The inset shows collapsed data of steady state C(r) 
for L = 64, 128, 256 and 512 as a function of r/L; the scaling 
function has a cusp with a ~ 0.25. 

Finally we note that the 2-point correlation functions 
in the steady state of the SP model exhibit a scaling form 
inr/L and have the same cusp exponents as in the coars- 
ening regime (with £ being replaced by L). For the EW 
surface, the scaling curve shown in the inset of Fig. |23| 
exhibits a cusp with a ~ 0.5. The corresponding curve 
for the KPZ surface, shown in the inset of Fig. ^ also 
exhibits a cusp, with a ~ 0.25. The fact that rric < 1 
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in these curves, as for those in the coarsening regime 
of the SP model, is indicative of the fact that the pure 
phase which are particle rich also have holes in them. In 
this respect, the pure phases differ from their CD model 
counterparts. 

We have seen above that the FDPO of the sliding par- 
ticles in the SP model is qualitatively of the same type as 
in the CD models for the underlying surfaces. We mea- 
sured the average overlap O — (siCTi) to get a quantitative 
estimate of the extent of correlation between the sliding 
particles (holes) and the valleys (hills) of the underlying 
surface. We found that it is nonzero as we expected, e.g. 
for the EW surface O ~ 0.26 and 0.39 corresponding to Si 
being defined within CD2 and CDS models. The overlap 
is greater in case of CDS model, since the domains are 
most often smaller than L/2 and this matches with the 
fact that particles clusters are also of size < L/2. On the 
other hand, domains in the CD2 model can be almost as 
big as L. For the KPZ surface, O ~ 0.26 corresponding 
to the overlap between particles and the coarse grained 
depth variables {s^j's of the CDS model. 



V. ROBUSTNESS OF FDPO 
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FIG. 25. S/C is plotted against kC, with R = 5 (triangles), 
R^l {•) and R = 0.2 (o) for the EW and KPZ surfaces at 
t — 400 X 2®. For clarity of display, we have multiplied the 
data for the KPZ surface by a factor of 2. 



A similar evaluation of C was also done for the KPZ 
surface, and is also shown in Fig. |2^ (upper curves). The 
observed slope of —1.25 implies that the cusp exponent 
remains a ~ 0.25 for all of them. 

We conclude that the variation of update rates affects 
the degree of ordering but not the asymptotic scaling 
properties, as indicated by the constancy of the cusp ex- 
ponent a. 



We did several numerical tests to check the robustness 
of the fluctuation dominated ordered state for the sliding 
particle (SP) problem. 

(i) We explored the effect of varying the ratio R = P2/P1, 
the relative rate at which the particles get updated as 
compared to the surface. 

(ii) We allowed the possibility of a small but finite rate 
(92 7^ 0) of the particles to hop uphill on a local r slope. 

(iii) We made the overall slope nonzero, in the case of 
the KPZ surface. 

We found that FDPO stays with (i) and (ii), while it 
is lost with (iii). 

For the EW surface, with R — 0.2 (i.e. the surface 
moving 5 times slower than the sliding particles), we 
found that Qi remains close to but slightly larger than 
0.18, the value for i? = 1. We checked the correlation 
function C(r, t) in the coarsening regime, and found that 
it has a cusp as a function of r/£ with the exponent 
a ~ 0.5. For R = 5 (i.e. the surface moving 5 times 
faster), we found Qi ~ 0.15. The latter value indicates 
a lesser degree of ordering and this is also mirrored in 
the 2-point function C{r,t): the collapse of the data as 
a function of r/C occurs beyond a time which is greater 
than that for i? = 1, i.e. the scaling regime sets in much 
later. Nevertheless, at large enough times, the cusp ex- 
ponent is unchanged (a ~ 0.5). Figure ^ (lower curves) 
shows the log- log plot oi S / C versus kC for the three 
rates i? = 5, 1 and 0.2. All of them have slopes —1.5, 
which indicate a ~ 0.5. 
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FIG. 26. 5/£ is plotted against kC, for t = 400 x 2^, with 
a finite uphill hopping rate for both EW and KPZ surfaces. 
The data for the KPZ surface is multiplied by a factor of 2 
for clarity of display. 

So far we have considered uphill hopping rate to be 
strictly zero, i.e. P2/92 — 00. By allowing for (72 7^ 0, i.e. 
allowing for upward motion of particles, we saw that the 
FDPO persists, so long as p2 > q2- In Fig. |6[ we show 
S{k)/C as a function of kC, at a large time t for EW and 
KPZ surfaces respectively, when the ratio = 5. We 

find that the slopes are —1.5 and —1.25 in the two cases 
indicating that the values of the cusp exponents are still 
a ~ 0.5 and a ~ 0.25 respectively, for the two surfaces. 

This points to the universality of the value a ~ 0.5 
(EW) and a ~ 0.25 (KPZ) over a range of models with 
different values of R, and also with respect to varying 
P2/q2- 
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FIG. 27. For a tilted KPZ surface, the curves for C{r,t) as 
a function of r all overlap at different times t = 400 x 2" (with 
n — 0, 6), indicating that there is no growing length scale 
£{t). Thus tilt removes FDPO. 

We also investigated the effect of having an overall tilt 
of the KPZ surface. This leads to an overall movement 
of the transverse surface fluctuations, which arc the ana- 
logues of kinematic waves in particle systems ■ In 
the presence of such a wave, the profile of hills and val- 
leys of the surface sweep across the system at finite speed, 
and the particles do not get enough time to cluster. Con- 
sequently the phenomenon of FDPO is completely de- 
stroyed. In Fig. ^ we show C{r,t) as a function of r 
(there is no scaling by C{t)) for several t. The curves are 
independent of t, in the absence of coarsening towards a 
phase ordered state. 



VI. CONCLUSION 



In this paper we have discussed the possibility of phase 
ordering of a sort which is dominated by strong fluctua- 
tions. In steady state, these fluctuations lead to vari- 
ations of the order parameter of order unity, but the 
system stays ordered in the sense that with probability 
one, a finite fraction of the system is occupied by a sin- 
gle phase. The value of this fraction fluctuates in time, 
leading to a broad probability distribution of the order 
parameter. 

We demonstrated these features in two types of mod- 
els having to do with surface fluctuations — the first, 
a coarse-grained depth (CD) model where we could es- 
tablish these properties analytically, and the second a 
model of sliding particles (SP) on the surface in ques- 
tion. For these models we found that besides (a) the 
broad probability distribution of the order parameter 
(which we may take to be the defining characteristic of 
FDPO), the steady state was also characterised by (b) 
power laws of cluster size distributions and (c) cusps in 
the scaled two-point correlation function, associated with 
the breakdown of the Porod law. The connection between 
(b) and (c) was elucidated using the independent interval 
approximation. Further, an extremal statistics argument 
showed that the largest cluster drawn from the power law 
distribution is of the order of the system size; this implies 



a macroscopic ordered region, so that within our models, 
properties (a) and (b) are connected. 

There are several open questions. Does fluctuation- 
dominated phase ordering occur in other, completely dif- 
ferent types of systems as well ? Are properties (b) and 
(c) necessarily concomitant with the defining property 
(a) of FDPO? Can one characterize quantitatively the 
dynamical behaviour in the FDPO steady state? 

Our model of particles sliding on a fluctuating sur- 
face relates to several physical systems of interest. First, 
it describes a new mechanism of large scale clustering 
in vibrated granular media, provided the vibrations are 
random both in space and time. Second, it describes a 
special case (the passive scalar limit) of a crystal driven 
through a dissipative medium, for instance a sediment- 
ing colloidal crystal |2|] . Finally, related models describe 
the formation of domains in growing binary films . It 
would be interesting to see if ideas related to FDPO play 
a role in any of these systems. 

It would also be interesting to examine fluctuating 
phase-ordered states in other nonequilibrium systems 
from the point of view of FDPO. For instance, in a study 
of jamming in the bus-route model studied in |2^, the 
largest empty stretch in front of a bus was found to be 
of order L, and it is argued that such a stretch survives 
for a time which is proportional to for a nonvanish- 
ing rate of arrival of the passengers. These features are 
reminescent of the behaviour of the CD and SP mod- 
els derived from the Edwards- Wilkinson model discussed 
above. However, more work is required to make a clear 
statement about FDPO in the bus-route model. 

In general, fluctuation-dominated phase ordering is ev- 
idently a possibililty that should be kept in mind when 
discussing new situations involving phase ordering in 
nonequilibrium systems, both in theory and in experi- 
ment. 

We acknowledge useful discussions with M.R. Evans, 
D. Dhar, G. Manoj, S. Ramaswamy and C. Dasgupta. 
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